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ABSTRACT 

The radio-quiet quasar SDSS J1536-I-0441 shows two broad-hne emission systems that Boroson & 
Lauer interpret as a candidate binary black- hole system with a separation of 0.1 pc (0.02 mas). From 
new VLA imaging at 8.5 GHz, two faint sources, separated by 0.97" (5.1 kpc), have been discovered 
within the quasar's optical localization region. Each radio source is unresolved, with a diameter of 
less than 0.37" (1.9 kpc). A double radio structure is seen in some other radio-quiet quasars, and the 
double may be energized here by the candidate 0.1-pc binary black- hole system. Alternatively, the 
radio emission may arise from a binary system of quasars with a projected separation of 5.1 kpc, and 
the two quasars may produce the two observed broad-line emission systems. Binary active galactic 
nuclei with a kpc scale separation are known from radio and X-ray observations, and a few such system 
are expected in the Boroson & Lauer sample based on the observed clustering of quasars down to the 
10 kpc scale. Future observations designed to distinguish between the 0.1 pc and 5 kpc scales for the 
binary system are suggested. 

Subject headings: quasars: individual (SDSS J1536-I-0441) — radio continuum: general — X-rays: 
general 



1. MOTIVATION 

The quasar SDSS J153636. 22+044127.0, at a redshift 
0.388, shows two broad- line er nission systems separate d 
in velocity by 3500 km s~^ (jBoroson fc Lauen l2009f ) . 
This unique quasar is interpreted as a binary black- 
hole system with a separation of 0.1 pc and masses of 
10^ '^ Mq and 10®'^ M©. The sub-parsec scale is signifi- 
cant, as it suggests that this close binary system shares a 
common narrow-line emission region and has also solved 
its so-called final parsec problem. For the assumed flat 
cosmology^, the 0.1-pc binary subtends 0.02 mas. An al- 
ternative interpretation for SDSS J153636. 22-1-044127.0 
(SDSS J1536-t-0441 hereafter) is that it represents two 
unrela ted quasars viewe d , by c hance, along similar sight 
lines. IBoroson &: Laueil ()2009l ) use the optical localiza- 
tion region - a circle of radius 1" - to rule out this alter- 
native interpretation with a probability of 0.0032. 

New high quality spectroscopy by IChornock et al.l 
(|2009D revealed a broad feature redshifted by 3500-4500 
km in Ha and H/3, which raises a third option that 
SDSS J1536-h0441 is a "do uble peaked emitter" (DPE) 
such as seen in Arpl02B jH alpern fc Filippenkol [Mi) 
and 3C 390.3 (jEracleous fc H alpern 1994). the pecuHar 
emission profile in DPE may originate from a thin Kep- 
lerian gaseous disk (also proposed for SDSS_J1536+0441 
by Ga skell [2009]). However, as noted bv lChornock et aTl 
(|2009f ). DPE show two broad components of comparable 
strength, unlike seen here. 

It is clearly important to observe this quasar with sub- 
arcsecond resolution, for two reasons. First, the 0.1- 
pc binary hypothesis can be ruled out if two widely- 
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separated sources are detected. Second, if observations 
with sub-arcsecond resolution detect only a single source, 
then the random-projection estimate would drop below 
0.0032, serving to strengthen the 0.1-pc binary hypothe- 
sis. 

SDSS J1536+0441 is not detected in 1.4-GHz sky sur- 
yeys, so its fiux den sity is less than 1 mJy at 5" resolution 
(IWhite et al.lll997D and l ess than 2.5 mJy at 45" resolu - 
tion (t Condon et al lll998 \ Following H^esic et all (|2002l ). 
combining the 1 m Jy up per li mit at 1.4 G Hz with the 
i-band magnitude ijBoroson fc Laueil |2009[ ) means that 
SDSS J1536-I-0441 has a radio-to-optical ratio less than 
unity, implying that it is a radio-quiet quasar. Fortu- 
nately, such a classification does not exclude it from be- 
ing detected at radio frequencies and observed with sub- 
arcsecond resoluti on. §[2]of this Letter reports new VLA^ 
([Thompson et al.lfl98Q ) imaging of SDSS J1536+0441 at 
sub-arcsecond resolution, leading to the discovery of two 
sources separated by 0.97" (5.1 kpc). The implications 
of this discovery are explored in § [S] § |4] closes with a 
summary of this work and suggestions for future direc- 
tions. 

2. VLA IMAGING 

The B configuration of the VLA was used, under pro- 
posal code AL738, to observe SDSS J1536-h0441 near 
transit on UT 2009 February 17 and 20. Observa- 
tions were phase-referenced to the calibrator J1534-I-0131 
whose positional accuracy was less than 2 mas. The 
switching angle was 3° and the switching time was 280 s. 
The center frequency was 8.4601 GHz, abbreviated as 
8.5 GHz hereafter. The a priori pointing position for the 
quasar was centered 3" north of the SDSS name-derived 
position to avoid any phase-center artifacts. Observa- 
tions were made assuming a coordinate equinox of 2000. 

Operated by the National Radio Astronomy Observatory, 
which is a facility of the National Science Foundation, operated 
under cooperative agreement by Associated Universities, Inc. 
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Data were acquired with a bandwidth of 100 MHz for 
each circular polarization. Observations of 3C 286 were 
used to set the amplitude scale to an accuracy of about 
3%. The net exposure times for SDSS J1536+0441 were 
6020 s and 6080 s on 2009 February 17 and 20, respec- 
tively. Twenty-five antennas provided data of acceptable 
quality. 

The data were calibrated using the 2009 Dec 31 re- 
lease of the NRAO AIPS software. No polarization cal- 
ibration or self-calibrations were performed. After cali- 
bration, each day's visibility data for SDSS J1536-I-0441 
were concatenated. The AIPS task imagr was applied 
to the concatenated data to form and deconvolve images 
of the Stokes / emission. A variety of weighting schemes 
were applied to the visibility data, converging on the im- 
age given in Figure 1 that optimizes the balance between 
sensitivity and resolution. 

Figure 1 shows two radio sources, labeled VLA-A and 
VLA-B. Elliptical-Gaussian fits to those sources in the 
image plane yielded the following integrated flux den- 
sities, positions, and 1-dimensional position errors: for 
VLA-A, S = 1.17±0.04mJy, a(J2000) 15''36™36!222, 
(5(J2000) = -|-04°41'27'.'06, and ctvla = O'.'l; for VLA- 
B, S ^ 0.27 ± 0.02 mJy, a( J2000) = 15''36™36!288, 
(5(J2000) = -f04°41'27'.'09, and ctvla = O'.'l. For each 
source, the flux-density error is the quadratic sum of the 
3% scale error and the fit residual, while the position 
error is the quadratic sum of a term due to the phase- 
calibrator position error (less than 0.002"), the signal- 
to-noise ratio (S/N) (less than 0.02"), and the phase- 
referencing strategies (estimated to be 0.1"). The im- 
age fits indicated that each source was unresolved and, 
given the high S/N data, this was taken to imply a di- 
ameter of less than 0.37", corresponding to half of the 
geometric-mean beam width at FWHM quoted in Fig- 
ure 1. Sources VLA-A and VLA-B have a summed flux 
density of 1.44 ±0.05 mJy. The image fits also yielded a 
relative source separation of 0.97" ± 0.03". 

3. IMPLICATIONS 

From the new VLA data, the 8.5-GHz emission from 
SDSS J1536+0441 consists of two sources, VLA-A and 
VLA-B, with radio luminosities, Lji = vLy, at 8.5 GHz 
of 5.2 X 10*° ergs s~-^ and 1.2 x 10*° ergs s~^, respectively. 
Each source is unresolved, with a diameter of less than 
0.37" (1.9 kpc) and, as Figure 1 sho ws, lies within the 
quasa r's optical localization region (|Boroson fc Laueil 
[2009t) . VLA-A and VLA-B are separated by 0.97" (5.1 
kpc). These could be two related radio sources, both 
energized by the candidate 0.1-pc binary system. Alter- 
natively, VLA-A and VLA-B could be independent radio 
sources originating from a binary quasar system, with 
a projected separation of 5.1 kpc. The implications of 
these two interpretations are explored separately below 
in §|3ll]and §[321 As mentioned in § [I1 SDSS J1536-h044 1 
is a radio-quiet quasar as defined bvl lvesic et al.l (|2002l ). 
and its properties will be analyzed within that context. 
Finally, in §|33]the available data for SDSS J1536-h0441 
will be used to explore a specific framework for radio - 
quiet quasars recently proposed bv lLaor fc Behail(|2f)0l . 

3.1. VLA-A and VLA-B Powered by a 0.1-pc Binary 

The radio emission from VLA-A and VLA-B could 
represent two related sources, with a projected separa- 



tion of 0.97" (5.1 kpc) and with both radio sources ul- 
timately being energized by the candidate 0.1-pc binary. 
When imaged at sub-arcsecond resolution, other radio- 
quiet quasars at similar redshifts and luminosities are 
known to exhibit double, triple, and linear radio struc- 
tures on scales of a few kiloparsecs (jKellermann et aP 
119941 : iKukula et al.l |1998[) . making a radio-double sce- 
nario plausible for J1536-I-0441. 

If either VLA-A or VLA-B coincide with the ener- 
gizing 0.1-pc binary, then the localization of that bi- 
nar y could be improved f rom the optical diameter of 
2" (jBoroson fc Lauer1l2009f ) to the radio di ameter of less 
than 0.37". However, by analogy with the iKukula et al.) 
(jl998) radio doubles with accurate optical astrometry, 
the 0.1-pc binary in J1536-I-0441 could be located be- 
tween VLA-A and VLA-B. Therefore a conservative 
value for the radio localization area is a region extend- 
ing about 1" east-west and less than 0.8" north-south 
(the major axis of the Gaussian restoring beam), lead- 
ing to a radio localization area of less than 0.8 arcsec^. 
The area of the optical localization region is 3.1 arcsec^ 
(jBoroson fc Lauerl I2009D . so the new 8.5-GHz imaging 
has improved the localization by a factor of at least 3.9. 

As mentioned in §[I1 SDSS J1536 +0441 was less tha n 
1 mJy at 1.4 GHz at 5" resolution (|White et al.lll997D . 
In contrast, the summed flux density for VLA-A and 
VLA-B is 1.44 ± 0.05 mJy at 8.5 GHz (§ [2]). Compar- 
ing these values suggests that the overall spectral index 
for SDSS J1536-I-0441 is flat, or even rising, with fre- 
quency. This could indicate that VLA-A, the stronger 
of the two 8.5-GHz sources, is sufficiently compact to be 
synchrotron self-absorbed and thus resemble components 
in some other radio-quiet quasars ( Baryainis et al.lfl996l : 
IKukula et allll998l : lUlvestad et"alll2005( l. Such compact 
emission could also hint that VLA-A marks the location 
of the candidate 0.1-pc binary, implying an even smaller 
radio localization for it and also strengthening its simi- 
larities to the radio galaxy 402-1-379 that hosts a 7-pc 
binary (|Rodriguez et al.|[2006l ). Some radio-quiet qua sars 
do exhibit time variability (e.g.. lBarvainis et al.ll2005[ ). so 
this inference of a flat or rising overall spectrum for VLA- 
A is weakened by not having simultaneous measurements 
at 1.4 and 8.5 GHz. But if the spectral index estimate 
is inexact due to variability on a decade timescale then, 
from causality arguments, the inference about small-scale 
radio emission still stands. 

3.2. VLA-A and VLA-B Powered by a 5-kpc Binary 

The detection of two point sources at radio frequen- 
cies also raises the possibility that the emission originates 
from two quasa rs 0.97" apart, or at a p rojected separa- 
tion of 5.1 kpc. iBoroson fc LaueJ (|2009( l noted that the 
probability for such a random projection in their sam- 
ple is 0.0032. Therefore, the two quasars are most likely 
not due to a random projection, but are likely physically 
related, i.e. a binary quasar system. 

A remarkable radio-loud case of a binary system of 
active galactic nuclei (AGNs) is seen in 3C 75, where 
two systems of two-sided jets emanate from two close 
point sources with a projected separation of ~ 7.5 kpc 
(jOwen et al.lll985f l. Compact binary AGNs were also 
revealed with Chandra observations of nearby systems. 
Particularly clear cases are NGC 6240 with a '--^ 1 kpc 
separation (|Komossa et al.l [2003D . and Mrk463 with a 
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~ 3.8 kpc separation (|Bianchi et al.ll2008D . Thus, SDSS 
J1536-I-0441 may be another example of a compact bi- 
nary AGN system. 

A systematic study of the abun dance of binary quasar s 
in the SDSS was carried out by iHennawi et aU (|2006l ). 
who found a projected correlation function of the form 
(-Rprop/0.43Mpch-i)-i-4« on scales of 10-40 kpc in 
proper length (proper length is used given the absence 
of a Hubble flow on these small scales). Extrapolating 
to 5 kpc (using h = 0.71) we get a correlation func- 
tion of 1175, i.e. an observed surface density 1175 larger 
than expected for random proje ction. Using the random - 
projection estimate of 0.0032 (jBoroson fc LaueJ I2009D . 
the expected observed number is actually 3.76, i.e. 3 to 
4 such binaries are expected. The study of Hennawi et 
al. is based on quasars with a mean z ~ 1.5 due to the 
smaller number of quasars at lower z, and it does not 
extend down to 5 kpc due to the fiber angular resolution 
limit. However, the study of Hennawi et al. indicates 
that the probab ility to find a 5 kp c binary quasar in the 
sample used by iBoroson fc Laueil ([2009) is of the order 
of unity. 

Each of the two quasars in the binary may be con- 
tributing its broad-line emission system to the total light, 
producing the double broad-line system discovered by 
[Boroson & Laucr (2009). Thus, rather than having a bi- 
nary black- hole system on a scale of 0.1-pc, we may have 
a binary quasar system on a 5-kpc scale, most likely re- 
siding within two strongly interacting galaxies. 

The velocity separ ation of 3500 km s^^ 
(HE roson fc Laueri 120091) is rather large, but not 
implausible in a c luster of galax i es. A bout half the 
clusters studied by iCarlberg et al.l (|1996D show such an 
extent of velocities. The maxim um velocity differences 
in the quasar binaries studied bv ^Henn awi et al"] ()2006f ) 
is 1870 km s~^, but that study imposed a cap of 2000 
km on the binary velocity separation. At a projected 
relative velocity of 3500 km s~^, two galaxies cannot 
form a bound system, so the term binary quasar here 
does not refer to a physically bound binary system. 

3.3. The Radio/X-ray Connection 

For some small X-ray and radio selected samples 
of radio-quiet AGNs, a link between the radio and 
X-ray emission has been suggested by the correla- 
tion between the rad i o and the X-ray luminosities 
(iBrinkmann et al.ll2000l : ISalvato et a"l]|2004t IWang et all 
I2006D. While intri guing, such a radio/X-ray connection 
should be verified by using an unbiased survey. The PG 
quasar sample is a com plete optically-selected sample 
(jBoroson fc GreenI 119921), making it i ndependent of ra- 
dio and X-ray biases. iLaor fc Behail (j2(308l ) have used 
these radio-quiet PG quasars to demonstrate that (a) 
the radio and the X-ray luminosities are correlated over 
a large range of AGN luminosity and (b) the correla- 
tion follows Ln/Lx ^ 10~^, the well- established cor- 
relation for coronally active cool stars (jGuedel fc Bend 
[l993), where Lr = vL^, at 5 GHz and Lx is in the 0.2- 
20 keV band. For cool stars, the L^J Lx ~ 10~^ relation 
is accepted as a manifestation of coronal heating by en- 
ergetic electrons following magnetic reconncction, that 
subsequently g ives rise to X-ray em ission. By analogy 
with cool stars. iLaor fc Behail ()2008f ) conjecture that the 
radio emission in radio-quiet AGNs may also be related 



to coronal magnetic activity. 

This coronal framework can be tested for SDSS 
J1536-I -0441 by making use of its r ecent Swift observa- 
tion bv lArzoumananian et al.l ()2009D . carried out in 2009 
Februar y 4 and 5, just two weeks befo re our VLA obser- 
vations. lArzoumananian et al.l (|2009[ ) measure a 0.5-10 
keV luminosity of 5 x lO^** erg s~^, with a spectral slope 
of -1.5. This gives z/L^ = 2.3 x 10"^ ^ erg s~^ at 1 keV, and 
thus Lx of 1.4 X 10**^ erg s""'^ (see ILaor fc Behml (|2008f ) 
for the conversion of vL^, at 1 keV to Lx)- The total 
radio luminosity we find at 8.5 GHz is 6.4 x 10*" erg s~^, 
which extrapolates to = 8.3 x 10'"' erg s~^, assum- 
ing a spectral slope of —0.5. We therefore get Lr/Lx = 
5.9 X 10~^ which falls within the range of ratios seen 
bv ILaor fc Beha^ ([2008) for the radio-quiet PG quasars. 
Thus, SDSS J1536-h0441 follows the radio/X-ray relation 
of optically selected radio quiet quasars. 

4. SUMMARY AND FUTURE DIRECTIONS 

The radio-quiet quasar SDSS J15 36-I-0441 has two 
broad -line emission systems that IBoroson fc Laueil 
(|2009f ) interpret as a candidate binary black-hole system 
with a separation of 0.1 pc (0.02 mas). Our new VLA 
imaging at 8.5 GHz reveals two sources, separated by 
0.97" (5.1 kpc), within the quasar's optical localization 
region. Each radio source has a diameter of less than 
0.37" (1.9 kpc). 

Other radio-quiet quasars do exhibit double structures, 
suggesting that the radio double in SDSS J1536-f0441 
could be energized by the candidate 0.1-pc binary. Al- 
ternatively, the radio emission may arise from a binary 
system of quasars with a projected separation of 5.1 kpc, 
and those two quasars may be responsible for the two ob- 
served broad-line emission systems. Binary AGNs with 
kpc-scale separations are known from radio and X-ray 
observations, and a few such system are expected in the 
IBoroson fc Laueil (|2009f ) sample based on the observed 
clustering of quasars down to a scale of 10 kpc. 

I nterestingly, the B aim er emission line profiles shown 
by 'Chornock et al. ( 2009) arc noticeably different than 
those in Boroson fc Lauer (2009). For H/3 the the peak 
of the r component is about 40% of the peak of b com- 
ponent, as measured above the "valley" betw e en th e r 
and b components, while in IBoroson fc Laueil (|2009l ) it 
is about 60%. Similarl y, for Ha the r peak i s about 
75% of the b peak in IChornock et all (|2009D . while 
in IBoroson fc Laueil (|2009D it is about 95%. These 
could be due to profile variability, but it could also 
be due to different contributions of two point sources 
to the total spectrum, due to differ ent slit posit i ons in 
the Boroson & Lauer (2009) and Ch ornock erall ()2009l ) 
spectra. 

If the 5-kpc binary interpretation is correct then 
ground-based adaptive optics imaging in the visible 
band, with sub-arcsecond resolution, should clearly show 
two point sources, and spectroscopy of those sources 
should exhibit two different spectra, each one with a sin- 
gle broad emission-line system. Deep optical imaging on 
arcsecond scales may also reveal the nature of the host 
galaxy, and whether strongly interacting galaxies are in- 
volved. 

Deep optical imaging on arcminutc scales can establish 
whether or not SDSS J 1536-I-0441 resides in a cluster. 
IBoroson fc Laue^ (|2009l) noted two nearby galaxies with 
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similar photometric redshifts, which may indicate the 
presence of such a cluster. However, a cluster is likely in 
both the 0.1-pc and 5-kpc binary interpretations, as a bi- 
nary black-hole also requires an earlier merger. The only 
difference between the two interpretation is the timescale 
since the merger occurred, which can be a fraction of the 
Hubble time for a 0.1-pc binary. On such a timescale the 
immediate environment may relax, but the larger scale 
cluster environment evolves independently of the merger 
process. 

Follow-up radio observations at other frequencies and 
epochs, as well as with higher angular resolutions, can 
help distinguish between the 0.1-pc and 5-kpc binary in- 
terpretations. Specifically, demonstrating that each ra- 
dio source has a flat spectrum, is time- variable, and re- 
mains compact at higher resolution would strongly sup- 
port the 5-kpc binary interpretation. In contrast, in the 
0.1-pc binary case, one (or both) of the radio sources 
would be expected to be steady in time, have a steep 



spectrum, and begin to show an outflow-like structure 
when examined with higher resolution; a flat-spectrum 
source paired with a steep-spectrum source would sug- 
gest a core-jet morphology. 

If the 5-kpc binary interpretation is correct, then SDSS 
J1536-I-0441 may form the first known compact binary 
quasar of two luminous broad-line systems, unlike the 
earlier X-ray-discovered compact systems where one or 
both of the AGNs are partly obscured, or broad emission 
lines are missing altogether. 
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Fig. 1. — VLA image of Stokes / emission from SDSS J1536+0441 at a frequency of 8.5 GHz and spanning 6" (31 kpc). SDSS position 
is marked with a circle of radius 1". Labels identify sources VLA-A and VLA-B. Uniform weighting with a robustness parameter of zero 
was used, giving an rms noise of 0.013 mjy beam~^ (1 a) and Gaussian beam dimensions at FWHM of 0.75" times 0.72" at a position 
angle of 1.3° (hatched ellipse). Geometric-mean beamwidth is 0.73" (3.8 kpc) at FWHM. Contours are at —6, —4, —2, 2, 4, 6, 8, 10, 
12, ... 20 times 1 a. Negative contours are dashed and positive ones are solid. Image peak is 1.15 mJy beam~^. Linear gray scale spans 
—0.05 mJy beam~^ to 1.0 mJy beam~^. 



